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The status of lattice calculations of heavy quark phenomenology is reviewed. Particular em- 
phasis is placed on the understanding and control of the calculational uncertainties. The ensuing 
implications for constraining the CKM matrix elements are discussed. 



I. INTRODUCTION 



X 



^\ , The calculation of hadronic matrix elements plays a vital role in determinations of the CKM matrix elements and 
ON ' overconstraining the unitarity triangle of the Standard Model. In particular, the CKM matrix elements that are 

combined to make up the sides of the triangle can be determined from neutral meson mixing and exclusive leptonic 

and semileptonic decay processes in the following way 

experimental measurement of 1 f known \ f nonperturbative 111/ I 2 (1) 

decay rate or meson mixing J 1 factors J 1 form factor or decay constant J ' CKM > ' <• > 

\ycKu\ can be determined by combining experimental measurements of exclusive decay rates with theoretical cal- 
culations of the nonperturbative contributions. These calculations can be done with a number of methods including 
' lattice QCD. Because the lattice formalism of QCD is model independent and systematically improvable it should 
C*~) ■ offer some of the most reliable results. Currently, the uncertainty in many of the theoretical and experimental inputs 
is large enough that unitarity is not really tested. Table | shows the decay modes to be discussed here and the asso- 
ciated CKM matrix elements. Indeed Vtd and V u b are amongst the least well determined components with much of 

ON ■ 

the error coming from the nonperturbative contributions (!]]. However, the next generation of B-physics experiments 
1 will reduce the experimental uncertainty considerably and to take full advantage of this the theoretical uncertainties 



must also be reduced. 

The form factors and decay constants of Equation [l] are extracted from hadronic matrix elements, which can 



i 

p : 

qj " be calculated directly in lattice QCD. Lattice calculations of heavy quark quantities began in the late '80s and it 
• was quickly realised that they could provide important inputs to phenomenologically interesting parameters of the 
^ , Standard Model. Over the years the field has matured enormously as a result of a deep understanding of the particular 
systematic errors involved and a concerted effort to control or remove them, as described in Section II. The largest 



remaining uncertainty in the numerical results is the quenched approximation. As this is removed in the next several 
years the lattice results will be the results of QCD and not a model. The aim of this paper is to give a progress report 
for heavy quark physics from lattice QCD emphasising the new results where attention has been paid to the particular 
systematic uncertainties involved. The state-of-the-art leptonic decay constants are described in some detail and a 
brief report on some new ideas for semileptonic physics is included. A more extensive review can be found in Ref. [Q . 

TABLE I. The CKM matrix elements and nonperturbative (lattice) calculations which are discussed in the text. 
CKM Matrix element Process Lattice Calculation 

I Vtd | B)i — Bj mixing f By /B d 

\V ts \ B° s - B s ° mixing f Bs 

|K*| B ► nil/ f±^(q 2 ) 

\V c b\ B -> Dlv h±(io) 

B~^D*lu h Al (oj) 



*Fermilab is operated by the Universities Research Association under contract with the US Department of Energy 
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II. PROGRESS IN HEAVY QUARK PHYSICS 



The evolution of lattice calculations using heavy quarks is shown very nicely in a plot made by Andreas Kronfeld 
for the Heavy Quark '98 Workshop g]. Figure [l] shows the lattice determinations of Jb from the earliest calculations 
to the most recently published ones, demonstrating how calculations of fg have matured over the years. Advances 
in computing power and our understanding of heavy quark physics mean the central value and error bars have 
settled down remarkably. The last four points are the most recent values of Jb from the JLQCD Q, Fermilab ||, 
NRQCD § and MILC @ collaborations respectively. I also note the UKQCD ||J|] and APE @ collaborations have 
new determinations of fs which are not included in this plot but which will be included in the text. 
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FIG. 1. fs vs time-ordered journal reference, as in Ref. M. Symbols represent the different methods for heavy quarks used: 
V is extrapolation from m c h — ■> mi,; □ is mq — > oo; Q i s interpolation [m c hi oo]; A is NRQCD and O is the Fermilab approach. 



The remainder of this section details the dominant uncertainties in a lattice calculation of heavy quark matrix 
elements and how they are treated in different approaches. 



A. O(amq) discretistion errors 

Present day computing resources mean lattice calculations are generally done in a box of size < 3 fm and with 
a lattice spacing a ~ 0.07 — 1.2 fm. Simulating realistically heavy quark masses at these lattice spacings must be 
approached carefully due to discretisation errors proportional to the heavy quark mass. A brute force reduction of a 
is extremely costly, requiring orders of magnitude more computing power so improving lattice actions is a necessary 
step. One approach is to use an improved light-quark action, with the improvement coefficients of the action and 
operators calculated nonperturbatively. Simulations are done with quark masses around that of charm. Heavy quark 
effective theory is then used to guide an extrapolation from this region to that of the bottom quark. In some cases 
results of a calculation in the static approximation can be used - allowing an interpolation in quark mass instead of 
an extrapolation. Both UKQCD and APE collaborations use this method. However, this extrapolation to a region 
where O(am) errors are not under control introduces a systematic error which is hard to quantify. 

A different approach is to consider the B system as nonrelativistic: for b quarks (v/c) 2 « (0.3GeV/5.0GeV) 2 . 
Relativistic momenta are excluded by introducing a finite cut-off such that p ~ itiqv <C mq. Then, \p\/niQ <C 1 and 
the QCD Lagrangian can be expanded in powers of I/tuq pH] . This has been extremely successful for B physics, 
where the quark mass is large, so the expansion is convergent. Results from the GLOK group using NRQCD are 
shown here. The theory is nonrenormalisablc which means a formal continuum limit (a — *■ 0) docs not exist so results 
must be obtained at finite lattice spacing. Lattice artefacts are removed by including higher orders in I/tuq and a. 
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The final approach is that of the Fermilab group |T^] . It is a re- interpretation of the relativistic action which identifies 
and correctly renormalises nonrelativistic operators present in the so called light-quark action. Discretisation errors are 
then O(clAqcd) and not O(amQ). Also the existence of a formal continuum limit means a continuum extrapolation 
is possible. The result is an action valid for arbitrary am and for which the systematic errors can be understood and 
controlled thereby allowing simulations at the b quark mass. This approach is used by Fermilab, JLQCD and MILC. 

B. Lattice spacing dependence 

If the large mass-dependent errors are under control the remaining lattice spacing dependence can be quantified 
by repeating the matrix element calculation at three or more lattice spacings and extrapolating the result to the 
continuum (a = 0) limit. This has been done by the collaborations using the Fermilab approach in their calculations 
of heavy-light decay constants. So far it has only been done by one group for the semileptonic case In general, 
only a mild dependence on lattice spacing is observed and reliable extrapolations can be performed. 

C. Operator matching 

Matrix elements determined on the lattice must be related to their continuum counterpart as in 

(P'(k)\0^\P(p)) contmuum = Z^g 2 ){P'{k)\0^\P{p)) latUce . (2) 

The renormalisation factor, Z^ is usually calculated in perturbation theory and has some associated uncertainty. 
UKQCD and APE use a nonperturbative determination of both the Z^ and the improvement coefficients in their 
actions, as described in which reduces considerably the matching uncertainty. 

D. Quenching 

Most lattice calculations are done in the quenched approximation (omitting light quark loops) as a computational 
expedient. The uncertainty introduced by this approximation is difficult to quantify. There are now however, some 
partially unquenched calculations of heavy quark quantities and these have been used by many collaborations to 
estimate the effect of quenching on their results. 

III. LEPTONIC DECAY CONSTANTS 

The calculation of heavy-light leptonic decay constants is an important application of lattice techniques for two 
reasons. Firstly, the decay constants themselves are important phenomenological parameters. In particular, the B 
meson decay constant, which has not been measured experimentally, is an input in determinations of the CKM matrix 
element Vta through — B^ mixing, 

f Bd B Bd \V td V t * h \ 2 . (3) 

Arrid has been measured experimentally to sa 4% |l| and the factor in square brackets is also well determined. 
Interestingly Vtd is known to only « 25% accuracy and the bulk of that uncertainty comes from the nonperturbative 
input f B Bs d - The bag parameter, BB d is also calculated on the lattice but will not be discussed here, some recent 
reviews are in Refs. [0 and S. 

Secondly, the calculation of leptonic decay constants is an important test-ground for lattice heavy quark techniques. 
Control of systematic errors here inspires confidence in the more complicated determinations of form factors and partial 
widths in the semileptonic case. 
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Naturally a calculation of the B meson decay constant can also include the decay constants of the B s meson and the 
D system. A number of experiments have measured fo s so the lattice result can be compared with the experimental 
numbers for at least one heavy- light decay constant. The agreement, shown in Figure ^ should inspire confidence in 
the determination of fg . 
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FIG. 2. Comparing the world-average lattice result for fn s with the experimental data flEHl 



In recent years there have been a number of new calculations of Jb using different approaches described in Section 
II, however a common feature is a comprehensive treatment of the systematic uncertainties involved. Table |n] lists 
results from the most recent calculations of the heavy-light decay constants and one can see that they are in rather 
good agreement now. World averages for all these decay constants can be found in the review by Draper 

Since both NRQCD and the Fermilab approach work well at the B meson mass it is interesting to compare the 
results. Figure || shows that for f\f~M (the amplitude of the matrix element) at I/tub there is good agreement . 
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FIG. 3. A comparison of Fermilab and GLOK results, q" re/ers S to the scale at which the operator matching is done and the 
Fermilab result is in the range considered by GLOK. Results from GLOK are at f3 = 6.0, Fermilab results are in the continuum 
limit. 



To illustrate the care taken in most calculations, the systematic error quoted by the Fermilab group was arrived at 
by considering the following: 



(i) Excited state contamination is less than statistics ; 
(Hi) Lattice spacing dependence is less than statistics ; 



(it) Finite volume is less than statistics 

(iv) Heavy quark tuning is less than statistics 
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(v) Perturbation theory = 5% ; 

So the systematic error in this calculation is almost entirely from the perturbation theory used to match the lattice and 
continuum operators, as described in Section II. The pertubative calculation used is a one-loop, mass-dependent result 
from Aoki et al |2(]] for which the leading error is taken to be 0{a 2 s ). The groups who studied the dependence on lattice 
spacing found it to be gentle, examples of the continuum extrapolations are shown in Figure ^. The quenching error 
has been estimated from work by the MILC collaboration who produced the first results for an unquenched fs using 
a relativistic action. They found an increase in fg of w 10% at a finite lattice spacing. Thus a (+16MeV) estimate 
from Fermilab. Other collaborations have used this work to make similar estimates (some taking an increased central 
value and symmetric error bars) . Preliminary results for an unquenched fs after extrapolation to the continuum limit 
suggest that the effect may be larger than 10% and closer to 25% |2lJ]. This is supported by the work of the NRQCD 
group, who have recently estimated this effect to be +25% p2fl . 

The ratios of decay constants are also calculated. In a ratio many of the systematic errors may cancel leading to a 
more precise result. This can be used to place bounds on |Vtd/Vts| and (assuming three generations) |Vtd| by using 
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FIG. 4. Continuum extrapolations: Fermilab extrapolates fs/~M to a = (left plot) and MILC extrapolates fs (right plot). 
The Fermilab action has a milder dependence on a as expected since the action is more highly improved. 



TABLE II. Recent calculations of /s and }b b //s from different collaborations. The two preliminary results from UKQCD 



were obtained in separate analyses. 



Group 




a- 1 


Statistical Error 


Systematic Error 


Quenching 


/bs//b 




[MeV] 


[GeV] 


[MeV] 


[MeV] 


[MeV] 




FNAL § 


164 


a -> 


+ 14 
-11 


±8 


+ 16 
- 


1.13 ± I 


MILC § 


159 


a -> 


±11 


+ 22 

- 9 


+ 21 

- 


1.11±2 + |±3 


jlqcdJB 


173 


a -> 


±4 


±9 


±9 




GLOK |f 


147 


1.92 


±11 


±11 


+ 8 
-12 


1.20±4+ o 


UKQCD §§ (prelim.) 


176, 161 


2.64, 2.58 


± I ±16 








APE 


179 


2.75 


±18 


+ 26 
- 9 




1.16±4 
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IV. SEMILEPTONIC B MESON DECAYS 



The lattice treatment of semileptonic decays of heavy-light mesons is not yet as mature as that of the leptonic 
decays which is clear from the absence of a complete estimate of systematic errors a la /b ■ The good news is that this 
is underway by a number of groups and should be available soon. In particular, calculations using NRQCD and the 
Fermilab approach herald a new era of precision in lattice calculations of form factors because the mass-dependent 
uncertainties are decoupled from other effects as a result of working at the b quark mass. Only B meson decays will 
be discussed here but D — ► -nlv and D — ► Klv are also being studied [p 23 1. 



A. B 
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To determine |V^(,| the required parameters are the form factors, f + '°(q 2 ) given by 



(P'{k)\V,\P(p)) = [p + k-q 



nip — m. 



r 



: /V), 



(5) 



where q 2 is the momentum transfer ie. g M = (p — fc) M . These decays are trickier than the leptonic case because 
the final state hadron introduces an extra kinematic parameter - the momentum of the final state particle. Lattice 
calculations work best at small three-momenta because momentum-dependent errors are under control ie. close to 
9max — i m P ~ mpi) ■ However, it is traditional to quote form factors at q 2 = since at this point experiments can 
measure a slope and intercept for the form factor, so a large model-dependent extrapolation in q 2 is performed. In 
contrast, for D — > ttIv the entire kinematic range can be covered, eg. Ref. [p3[ . Figure || shows preliminary data from 
Fermilab for the form factors as a function of q 2 . For the first time the lattice spacing dependence of these matrix 
elements has been studied and for 300MeV < p n < 800MeV this is found to be mild (b|] . Finally I note that | V u b I can 
be determined without recourse to a (^-extrapolation by calculating partial widths from the differential decay rate 
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300MeV<p^<800MeV 

in a region where theory and experiment have reliable data, as in Figure [|. 
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FIG. 5. The q 2 dependence of the form factors by the Fermilab group for strange light quarks and a — > 0, compared with 
results from UKQCD and Light cone sum rules jjq ]. The lattice data only include statistical errors. The right-hand plot 
is the decay rate (less the momentun independent prefactors) versus pion momentum. The dotted lines define a region in pion 
momentum where it is believed the lattice calculation is most reliable and for which there are experimental data. 
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B. B — > Dlu 

The final topic is a new development in the calculation of form factors at zero recoil for the decay B — > D^*Hv. 
The matrix element is parameterised by two form factors, ft/ + '~'(l) which are combined in a physical form factor, 

Tb^d(1) = h B ^ D {\) - h B ^ D (l)(m B - m D )/(m B + m D ) . 
\V c b\ is determined from experimental data for ,F(l)|T4f,|. Although the shape of the form factor has been calculated 
successfully (see eg. Refs. 26,3]) lattice calculations of the absolute value at zero recoil were frought with difficulties 



both statistical and systematic. But it has recently been shown that by extracting the form factors from a ratio of 
matrix elements at zero recoil, the bulk of the uncertainties cancel leaving an extremely precise result for the form 
factor |27]] and therefore \V c b\- For h + (l) one constructs the ratio 

(D\V \B){B\V \D) h B - D {l)h^ B (\) 

(D\V Q \D){B\V \B) h B ^ D {l)h B - B {\) 1 + [ >l ■ 1 1 

There is a similar expression using double ratios for h~(l). Preliminary results [f2~7| show this is indeed very successful, 

= 1.069 ± 0.008 ± 0.002 ± 0.025. (8) 

In fact this method has the added advantage that at zero recoil it is a deviation of .F(l) from 1 that is measured so 
that effects of quenching and lattice spacing dependence which do not cancel in the ratio are expected to be small. 
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